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Bilayer graphene has attracted considerable interest due to the important role played by many- 
body effects, particularly at low energies. Here we report local compressibility measurements of a 
suspended graphene bilayer. We find that the energy gaps at filling factors v — ±4 do not vanish 
at low fields, but instead merge into an incompressible region near the charge neutrality point at 
zero electric and magnetic field. These results indicate the existence of a zero-field ordered state 
and are consistent with the formation of either an anomalous quantum Hall state or a nematic 
phase with broken rotational symmetry. At higher fields, we measure the intrinsic energy gaps of 
broken-symmetry states at v = 0, ±1 and ±2, and find that they scale linearly with magnetic field, 
yet another manifestation of the strong Coulomb interactions in bilayer graphene. 

PACS numbers: 73.22.Pr, 73.43.-f 



The low-energy dispersion of bilayer graphene can be 
described to first order by parabolic valence and con- 
duction bands that meet at the Fermi energy [1]. The 
charge excitations of this band structure are massive chi- 
ral fermions whose quantum Hall signature is different 
from that of both monolayer graphene [2, 3] and con- 
ventional two-dimensional electron gases. In the absence 
of interactions, the Landau level (LL) energy spectrum 
of bilayer graphene is given by En = ±huj c ^N(N — 1), 
where N is the orbital index, ft = h/2n (h is Planck's 
constant), uo c = eBjm* is the cyclotron frequency, e is 
the electronic charge, B is the magnetic field, and m* 
is the effective mass [4]. Plateaus occur in Hall conduc- 
tivity at a xy = 4Me 2 /h, where the factor of four is due 
to spin and valley degeneracy and M is a nonzero inte- 
ger, reflecting the eightfold-degenerate LL at zero energy, 
which is comprised of TV = and 1 orbital states. 

When magnetic field is large enough or disorder is suf- 
ficiently low, interaction effects such as quantum Hall 
ferromagnetism (QHF) [5] or magnetic catalysis [6] are 
predicted to open energy gaps and give rise to additional 
plateaus in the quantum Hall spectrum at intermediate 
filling factors. Recent transport measurements have in- 
deed revealed signatures of many-body effects in bilayer 
graphene [7-10] and the dependence of resistance on tem- 
perature and magnetic field was used to determine the 
magnitude of the interaction-driven energy gaps [7, 8]. 
More recently, it has been theoretically predicted that 
spontaneously broken symmetries will occur in bilayer 
graphene at zero magnetic field. The nature of the zero- 
field interacting phase is still under intense theoretical de- 
bate, with suggestions of spontaneous transfer of charge 
between layers [11, 12] ferroelectric domains [13], nematic 
ordering [14, 15], or the formation of an anomalous Hall 
insulator [11, 12, 16-19]. 

One way to distinguish between the various interacting 
states is to determine whether an energy gap is present at 
B = 0. Transport measurements can provide an indica- 



tion of gap size, but it is known that disorder can decrease 
the apparent transport gap relative to the true intrinsic 
gap [20, 21]. Therefore, it is desirable to directly probe 
electronic properties in a thermodynamic measurement. 
Here, we report the use of a scanning single-electron tran- 
sistor (SET) [22-25] to measure the local compressibility 
of a suspended bilayer graphene flake. Investigations into 
the compressibility of bilayer graphene were recently re- 
ported on unsuspended samples [26, 27] with particular 
attention paid to the gap induced by an electric field, but 
disorder in these systems was too large to observe the 
broken-symmetry states discussed above. Our measure- 
ment combines the high sensitivity afforded by an SET 
with the low disorder of suspended devices, allowing us 
to study electronic states that arise from Coulomb inter- 
actions and revealing the existence of an ordered state at 
zero field. A schematic illustration of the measurement 
system is shown in Fig. 1(a); for a full description of the 
measurement technique, see refs. [22, 23, 28]. The SET is 
capable of measuring changes in local electrostatic poten- 
tial with fiV sensitivity. As carrier density n is varied, 
changes in <!> directly reflect the changes in local chemical 
potential /i of the bilayer flake, so the scanning SET tip 
can be used as a local probe of inverse compressibility 
(incompressibility) d/a/dn. 

In order to establish our measurement technique, we 
first describe the behavior of our sample in the high-field 
regime. Figure 1(b) shows a typical measurement of in- 
verse compressibility taken at B = 2 T. Peaks in d/n/dn, 
caused by the low density of states between neighboring 
LLs, are apparent at v = 0, ±2 and ±4. The observation 
of incompressible regions at v = and ±2 is indicative 
of broken symmetries in the zero-energy LL, consistent 
with recent transport measurements [7-10]. The widths 
of the incompressible regions provide a measure of disor- 
der [25], and the full width at half maximum at v = 4 
is on the order of 10 10 cm -2 . This is more than ten 
times smaller than in unsuspended devices [29], and is 
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FIG. 1. (Color online) (a) Schematic illustration of the mea- 
surement setup, (b) Inverse compressibility as a function of 
carrier density at B — 2 T. Incompressible peaks at v — 0, ±2, 
and ±4 occur due to the decreased density of states between 
Landau levels, (c) Chemical potential as a function of carrier 
density at B = 2 T, obtained by integrating the data in b. 
Steps in chemical potential occur at v = 0, ±2, and d=4 as 
electrons begin to occupy the next Landau level. In panels b 
and c, experimental data are shown in blue and fits, based on 
Lorentzians at each filling factor, are shown in red. 



FIG. 2. (Color online) (a) Gap size, extracted from a 
Lorentzian fit, as a function of magnetic field for v — (blue), 
±1 (magenta), ±2 (green), ±4 (red), and ±8 (black). Squares 
correspond to electrons, and circles to holes. Gap size at all 
filling factors is well described by linear scaling with magnetic 
field, with fits given by the solid lines, (b) Incompressible 
peak width as a function of magnetic field for v — 0, ±2, ±4, 
and ±8 (same colors as in a). The v — peak is significantly 
narrower than the others. Peak width does not strongly de- 
pend on magnetic field. 



consistent with estimates from transport measurements 
of similar suspended bilayers [7]. 

Chemical potential as a function of carrier density, 
shown in Fig. 1(c) for B = 2 T, is obtained either by 
direct measurement or by integrating curves similar to 
that shown in Fig. 1(b). The steps in fi(n) at each 
filling factor provide a measure of the energy gaps A u 
between neighboring LLs. We have measured gap size at 
v = 0, ±1, ±2, ±4 and ±8 as a function of magnetic field, 
and the resulting data are shown in Fig. 2(a). Con- 
sistent with the expected behavior of En^A^ and Ag 
scale linearly with magnetic field, with magnitudes of 3.9 
meV/T and 2.8 meV/T, respectively. This linear depen- 
dence confirms that the range of carrier densities probed 
in our experiment lies within the energy regime where 
the bilayer graphene band structure is well-approximated 
by parabolic bands. It should be contrasted with cy- 
clotron resonance and compressibility studies performed 
at higher densities where the hyperbolic nature of the 
dispersion was apparent from the sublinear dependence 
of gap size on magnetic field [26, 30]. 

It is apparent from Fig. 2(a) that the energy gaps 
Aq and A2 corresponding to broken-symmetry states at 
v — and ±2 also increase proportionally with B, and 
have gap sizes of 1.7 meV/T and 1.2 meV/T, respec- 
tively. For v = 0, the data show good agreement with 
the activation energy E a = A /2 = 0.3 - 0.9 meV/T 



measured in transport experiments [7], as well as recent 
theoretical predictions [13, 31, 32]. For v = 2, activa- 
tion experiments on unsuspended flakes [8] yielded sig- 
nificantly smaller gap sizes and suggested a B 1 ! 2 depen- 
dence for the gap, which seems to conflict with our mea- 
surement. However, the error bars from the transport 
measurement are large enough that a linear fit to the 
data which passes through the origin is not inconceiv- 
able. Linear scaling with B for interaction-driven LLs 
is reasonable if one considers screening from higher or- 
bital LLs, whose energy separation is much smaller than 
the Coulomb energy for all experimentally relevant fields 
[13]. The linear scaling of Aq and A2 can therefore be 
understood as a result of the very strong Coulomb inter- 
actions in bilayer graphene. Finally, our sample exhibits 
energy gaps of less than 1 meV at v = zbl in high mag- 
netic fields. Linear scaling with a slope of approximately 
0.1 meV/T provides a reasonable fit to these data, but 
the gap sizes are too small to conclusively rule out B 1 ! 2 
dependence. 

We can use the measured gap sizes to determine the 
effective mass of bilayer graphene. Summing the mea- 
sured energy gaps for \v\ < 4 at each magnetic field and 
comparing with our expression for En, we obtain an es- 
timate that m* = (0.027 ± .002)ra e , in good agreement 
with Shubnikov de Haas measurements [33] . From a sim- 
ilar analysis using the gap sizes for |z/| < 8, we extract 
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an effective mass of (0.032 ± .002)m e , which is somewhat 
larger than the above value. The origin of this discrep- 
ancy is unclear, but it may be an artifact of our analysis 
procedure [28] . It is also important to note that for sam- 
ples with a single gate, density and electric field cannot 
be controlled independently, so all gaps at nonzero filling 
factors are measured in an electric field. Changes in the 
expected effective mass could result from deviations from 
the parabolic band structure or from an electric field ap- 
plied perpendicular to the flake [33, 34]. 

Figure 2(b) shows the widths of the incompressible re- 
gions at each filling factor as a function of magnetic field. 
The peak at v = is significantly narrower than the oth- 
ers [see also Fig. 1(b)], which may indicate that higher 
filling factors are subject to additional sources of disorder 
that do not affect the v — state. One possible explana- 
tion for this finding is the existence of variations in the 
effective magnetic field, which can be caused by ripples 
or strain [35-37]. 

We now discuss the behavior of our sample at small 
magnetic fields. Figures 3(a) and 3(b) show inverse com- 
pressibility as a function of density and magnetic field. 
We observe distinct incompressible peaks corresponding 
to quantum Hall states at v — ±4 that extend all the way 
down to B = 0, where they merge into an incompress- 
ible region at the charge neutrality point [Figs. 3(a)-(c)]. 
Surprisingly, the v = ±4 gaps do not vanish at low fields; 
in fact, they increase with decreasing field below 0.2 T, as 
shown in Fig. 3(d). The zero-field incompressible peak 
has a characteristic width of approximately 10 10 cm -2 
[Fig. 3(c)], and integration yields a step in chemical po- 
tential of nearly 2 meV. This increase in chemical po- 
tential is too large to be explained by a disorder-induced 
electric field, which would only lead [33] to gaps on the 
order of 0.1 meV for charge inomogeneity of 10 10 cm -2 . 
We observe qualitatively similar behavior at all positions 
along the flake [28]. It is worthwhile to note that our 
measurements give no indication of the negative incom- 
pressibility [38] at zero field, but are consistent with more 
recent predictions [39]. 

To further elucidate the origin of the low-field incom- 
pressible behavior, we have performed transport mea- 
surements on the same flake. Two-terminal resistance 
as a function of carrier density and magnetic field is 
shown in Fig. 4(a). Despite the relatively large jump 
in chemical potential that we observe in compressibility 
measurements, the resistance near the charge neutrality 
point at zero magnetic field is only a few k^. For com- 
parison, we note that the gap size of the v = state at 
B — 0.5 T is approximately 0.5 meV, but even for this 
small gap, a marked decrease in conductance is already 
apparent [Fig. 4(c)]. Therefore, the formation of an en- 
ergy gap with conduction mediated by activation is an 
insufficient explanation for the zero-field behavior that 
we observe. The derivative of conductance with respect 
to carrier density is plotted in Fig. 4(b). Several sharp 
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FIG. 3. (Color online) (a) Inverse compressibility as a func- 
tion of carrier density and magnetic field. Incompressible 
peaks occur at v — and ±4. Below about 0.2 T, peak height 
along v — dz4 increases with decreasing field, culminating in 
an incompressible peak at zero field, (b), Zoom-in on the low- 
field behavior of the sample, taken at a different location on 
the flake from that in a. Distinct peaks along v — ±4 persist 
all the way to zero field, (c) Line cuts of incompressibility 
along the red arrows shown in panel b for magnetic fields be- 
tween B = and 0.175 T in steps of 0.035 T. Curves are 
offset for clarity. Data are shown in blue, and the red curves 
are two-Lorentzian fits, except for the zero-field fit, which is 
composed of only one Lorentzian. (d) Gap size at \u\ = 4 for 
electrons (green), holes (red), and their sum (cyan). Data for 
\B\ > 0.03 T is based on a two-Lorentzian fit. The blue cir- 
cles at low field describe the jump in chemical potential across 
the charge neutrality point, as modeled by a single Lorentzian 
fit. Solid black lines correspond to A4 = 3.9(£?[T]) meV, as 
derived from the data at high magnetic fields. 



lines with the same slope as the v — 4 state are apparent, 
and these conductance fluctuations are caused by local- 
ized states in the bilayer [40]. The localized states, which 
indicate the presence of an energy gap, persist all the way 
to zero field. 

The incompressible behavior and transport character- 
istics at B = indicate the presence of an interacting 
state. In the limit of parabolic bands, the non-vanishing 
density of states at the charge neutrality point means 
that even infinitesimally small electron-electron interac- 
tions can lead to correlated states [11-19]. The existence 
of a zero-field incompressible peak and the fact that gap 
size along v — ±4 does not vanish at low fields are con- 
sistent with two proposed interacting states, which we 
discuss below. 

One proposal that is consistent with our measurements 
is the formation of an anomalous Hall insulator at low 
electric and magnetic field [11, 12, 16-19]. In such a state, 
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FIG. 4. (Color online) (a) Two-terminal conductance as a 
function of carrier density and magnetic field. Quantum Hall 
plateaus at v — ±4 emerge at very low field, and the onset 
of the highly resistive state at v — is also apparent, (b) 
Derivative of conductance with respect to carrier density. Lo- 
cal peaks and minima with the same slope as v — ±4 are 
visible, indicative of localized states. These localized states 
persist down to zero magnetic field, (c) Conductivity along 
filling factors v — —4 (red), (blue), and 4 (green) as a func- 
tion of magnetic field. 



and because of the gapless nature of the spectrum, the 
conductance should remain relatively high, also consis- 
tent with our observations. 
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time-reversal symmetry is spontaneously broken and do- 
mains form where the flake is at either v = 4 or v — —4. 
It is reasonable that we observe incompressible behavior 
at both filling factors if we assume that the SET is too 
large to resolve individual domains. In this scenario, the 
B = conductance is dominated by edge state transport 
and should therefore remain at a few kQ per square, con- 
sistent with our findings. The fact that gap size along 
v = ±4 increases with decreasing field below about 0.2 
T [Figs. 3(a), 3(b) and 3(d)] can be understood to arise 
from the competition between the anomalous Hall phase 
and the standard quantum Hall gap as magnetic field is 
increased [18, 19]. 

So far, we have assumed parabolic bands and neglected 
trigonal warping. Trigonal warping modifies the band 
structure of bilayer graphene so that four Dirac cones 
emerge at low energies, leading to a 16— fold degenerate 
LL at zero energy for low magnetic fields. Effects of trig- 
onal warping should be apparent [41] at densities as high 
as 10 11 cm -2 , which is well within the experimentally 
accessible regime due to the low disorder in our sample. 
Our results, however, cannot be explained by trigonal 
warping in the single-particle picture because the incom- 
pressible peaks at low magnetic field in Fig. 4(a) follow 
v = ±4 rather than v = ±8. Recently, however, it was 
theoretically predicted [14, 15] that electron-electron in- 
teractions can break rotational symmetry and modify the 
dispersion into a nematic phase that is characterized by 
two Dirac cones. In this scenario, one would expect in- 
compressible peaks along v = ±4, with energy gaps that 
scale as ~ B 1 ' 2 . This means that gap size would re- 
main relatively large at low fields, as we observe. In this 
picture, the zero-field incompressible peak that we mea- 
sure can be ascribed to the vanishing density of states, 
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FIG. 1. Inverse compressibility as a function of carrier density at zero magnetic field. A sharp peak near the charge neutrality 
point is apparent, consistent with the measurements presented in the main portion of the paper. There is little variation in 
compressibility outside the central peak. 



METHODS 



Sample fabrication was performed according to the method described in ref. [1]. Graphene was deposited onto a 
Si/Si02 wafer by mechanical exfoliation, and suitable bilayer flakes were contacted using an electron-beam lithography 
step followed by thermal evaporation of Cr/Au contacts. Samples were then suspended using 5:1 buffered oxide etch, 
dried in a critical point dryer, and quickly transferred to the measurement system, an ultra-high vacuum scanning 
probe system with a 3 He refrigerator. All measurements were performed at base temperature of the system, 450 mK. 
Back gate voltage was limited to ±15 V to avoid collapse of the device. 

Scanning probe tips were fabricated using a fiber puller to generate a conical taper at the end of a quartz rod. Al 
(19 nm) was thermally evaporated onto either side of the quartz rod to generate the two leads of the single electron 
transistor (SET). Following a brief oxidation step, another 19 nm of Al was evaporated onto the end of the conical tip 
to serve as the island of the SET. For the measurements described in this paper, tip size was approximately 150 — 200 
nm and the tip was held between 50 nm and 200 nm above the flake. 

Compressibility measurements were performed using both AC and DC methods. In the AC scheme, a small AC 
bias is applied to the back gate in order to weakly modulate the carrier density (by about 1.5 • 10 9 cm -2 ) in the 
bilayer flake. The corresponding changes in SET current are measured and converted to an electrostatic potential, 
with the conversion factor obtained by simultaneously monitoring the SET response to a 1 mV AC potential applied 
directly to the sample. In the DC scheme, the current through the SET is maintained at a fixed value using a feedback 
mechanism that modifies the potential of the sample. The changes in sample potential provide a measure of the change 
in electrostatic potential caused by a given change in carrier density. In both cases, dji/dn and fi(n) are obtained 
because the electrostatic potential is related to the chemical potential of the flake by e<l> = \i. 

Due to the finite sample size, the SET is also affected in both schemes by fringing fields from the back gate [2, 3], 
which give rise to a constant offset in d/n/dn. We account for this parasitic capacitance by subtracting a constant value 
so that djii/dn = within each Landau level. This is necessarily an overestimate, because it corresponds to a delta- 
function representation of the Landau level, whereas our Landau levels are broadened by disorder. Therefore, our gap 
sizes should be taken as lower limits. Incompressible peaks at each filling factor are fit with a single Lorentzian, with 
amplitude and width as fitting parameters. The uncertainty in effective mass is calculated by adding in quadrature 
the uncertainty at each filling factor, as obtained by a linear regression. 
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FIG. 2. (a) Inverse compressibility as a function of carrier density and position along the flake at 0.1 T. The incompressible peaks 
at v — ±4 (centered at approximately ±10 10 cm -2 ) are consistently apparent, independent of position. The incompressible 
peak on the hole side is consistently stronger than that on the electron side. Some variations in the density at which the peaks 
occur is visible, providing an estimate of disorder, (b) Inverse compressibility as a function of carrier density and position along 
the flake at 0.2 T. The incompressible peaks at v — ±4 (centered at approximately ±2 • 10 10 cm -2 ) are again apparent along 
the entire flake. The electron-hole asymmetry is also apparent at this field, and the v — peak is just beginning to emerge. 

ZERO-FIELD COMPRESSIBILITY AT LARGE DENSITIES 

Figure SI shows the inverse compressibility of our sample at zero magnetic field, taken over a larger density range 
and at a slightly different position along the flake. Similar to the measurements presented in the main portion of 
the paper, a sharp incompressible peak is apparent at the charge neutrality point. As electron density is increased 
beyond this central peak, little change in compressibility occurs over the entire measurement range. For holes, inverse 
compressibility increases slightly as density is increased. The origin of this behavior is unclear, and the change is an 
order of magnitude smaller than the negative compressibility predicted by Kusminskiy et al [4] . It also appears that 
the compressibility is somewhat larger for holes than for electrons. Such behavior is consistent with the inclusion of 
the 74 hopping parameter, which increases the density of states for holes and decreases it for electrons [5]. However, 
a non-negligible 74 would also be expected to cause an asymmetry between the energies of electron and hole Landau 
levels, which we do not observe. 

SPATIAL SCANS AT LOW MAGNETIC FIELD 

It is apparent in Fig. 3 that the incompressible peak height as well as the gap size are larger for holes than for 
electrons along \v\ = 4 at low magnetic fields. Figures S2(a) and S2(b) show compressibility measurements taken as a 
function of density and position along the flake at B = 0.1 T and 0.2 T, respectively. These measurements show that 
the electron- hole asymmetry is consistent across the entire width of the sample. It is also worthwhile to note that the 
fluctuations in density at which the incompressible peaks occur can provide another measure of the disorder in our 
samples. Consistent with estimates made by other means, it appears that the charge inhomogeneity is less than 10 10 

— 2 

cm . 



ADDITIONAL MEASURES OF v = ±4 GAP MAGNITUDE 

Figs. S3(a) and S3(b) show slanted scans of inverse compressibility which follow a fixed carrier density range 
surrounding filling factors v = — 4 and 4, respectively, at magnetic field between and 1 T. At low fields, the second 
\u\ = 4 peak and the v = peak are visible as well. These measurements were performed at slightly different positions 
and heights compared to those presented in the paper. Figs. S3(c) and S3(d) show the extracted gap size at each 



4 





FIG. 3. (a) Inverse compressibility taken over a constant density range centered around v — —4 as a function magnetic field. 
Data was taken at a slightly different position and height compared to that presented in the main paper. At low fields, the 
v — 4 peak is visible as well, (b) Inverse compressibility taken over a constant density range centered around v — 4 as a function 
magnetic field. Data was taken at a slightly different position and height compared to that presented in the main paper. At 
low fields, the v — and —4 peaks are visible as well, (c) Gap size along v — —4 as a function of magnetic field, as fit by a 
single Lorentzian. For B < 0.1 T, the gap size may be an overestimate due the effect of the v — 4 incompressible peak, (d) 
Gap size along v — —4 as a function of magnetic field, as fit by a single Lorentzian. For B < 0.1 T, the gap size may be an 
overestimate due the effect of the v — 4 incompressible peak. 



respective filling factor as a function of magnetic field. We again observe that the gap size along v = ±4 does not 
vanish at low field, but increases with decreasing field below 0.2 T. For B > 0.2 T, gap size increases linearly with 
field, as expected. The qualitative agreement between this dataset and that presented in the main paper provides a 
second indication that our findings are not strongly dependent on position. 
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